LTH 862, DESY 10-007, arXiv : 1 1 . 3554 [hep-ph] 



PROCEEDINGS 

OF SCIENCE 



January 2010 



Higher-order predictions for splitting functions and 
coefficient functions from physical evolution kernels 



o 

C 

Hi 

o 

(N 



a,; 



> 

in 
cn 

o 
o 



X 



A. Vogt* 

Department of Mathematical Sciences, University of Liverpool, UK 



E-mail: Andreas . Vogt@liv. ac . uk 



S. Moch 

Deutsches Elektronensynchrotron DESY, Zeuthen, Germany 
E-mail: 



sven-olaf .moch@desy.de 



G. Soar 

Department of Mathematical Sciences, University of Liverpool, UK 
E-mail: |G.N.Soar@liv.ac.uk 



J. A.M. Vermaseren 

NIKHEF Theory Group, Amsterdam, The Netherlands 



E-mail: t68@nikhef.nl 



We have studied the physical evolution kernels for nine non-singlet observables in deep-inelastic 
scattering (DIS), semi-inclusive e + e~ annihilation and the Drell-Yan (DY) process, and for the 
flavour-singlet case of the photon- and heavy-top Higgs-exchange structure functions (F2, Fa) in 
DIS. All known contributions to these kernels show an only single-logarithmic large-x enhance- 
ment at all powers of (1 — jc). Conjecturing that this behaviour persists to (all) higher orders, 
we have predicted the highest three (DY: two) double logarithms of the higher-order non-singlet 
coefficient functions and of the four-loop singlet splitting functions. The coefficient-function pre- 
dictions can be written as exponentiations of 1 /rV-suppressed contributions in Mellin-rV space 
which, however, are less predictive than the well-known exponentiation of the ln k N terms. 



RADCOR 2009 - 9th International Symposium on Radiative Corrections (Applications of Quantum Field 

Theory to Phenomenology) 

October 25-30 2009, Ascona, Switzerland 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Higher-order predictions from physical evolution kernels 



A. Vogt 



1. Introduction: hard lepton-hadron processes in perturbative QCD 

We are interested in the structure functions in deep-inelastic scattering (DIS), the corresponding 
fragmentation functions in semi-inclusive e + e~ annihilation (SIA), and the cross section —da/dM^ 
for Drell-Yan (DY) lepton-pair production in hadron-hadron collisions (see Ref. [1] for a detailed 
introduction). These one-scale observables, here denoted by F a (x,Q 2 ), are generically given by 

F a (x,Q 2 ) = [C fliiW (0fe(Ai 2 ),M 2 /G 2 )®yi*(M 2 ){®//(M 2 )}]W + ^(l/G (2) ) • (1-1) 

Here Q 2 denotes the physical hard scale (e.g., Q 2 = Mfi for the DY case), and x the corresponding 
scaling variable. jJ. represents the MS renormalization and factorization scale (there is no need to 
keep them different here), and ® stands for the Mellin convolution. The parts of Eq. Ql.l[ ) in curly 
brackets only apply to the DY case, and summation over i {and j} is understood. 

At jj, 2 = Q 2 the expansion of the coefficient functions C a in powers of the strong coupling a s is 

C a ,i{x, 0£s) = (1 - 8 aL ) 8 iq 8{\-x) + a 5 c ( a l J (x) + a 2 ^ 2 ) (x) + a s 3 cg } (*) + ... (1.2) 

with a s = GC S / (4tt). As indicated by the first term of the r.h.s., of all cases we consider here only the 
longitudinal coefficient functions in DIS and SIA vanish at order a s °. The (spacelike) parton and 
(timelike) fragmentation distributions ff 1 are, of course, non-perturbative quantities. However their 
scale dependence is calculable perturbatively via the renormalization-group evolution equations 

d 



dlnju. 



f^ii 2 ) = (as(li 2 ))®f k (ji 2 )](x) , P(x,a 5 ) = Zi=o4 +1 P {l) (x) ■ d-3) 



Except for Fl in DIS and SIA, the terms up to ci'^ (x) and P (") (x) in Eqs. (1.2) and (1.3) define the 



N"LO approximations to Eqs. (|P|). Precise predictions including a sound numerical uncertainty 
estimate require, at least, calculations at the next-to-next-to-leading order (NNLO = N 2 LO). The 
same order is usually required for deducing structural features such as the ones discussed below. 

The NNLO coefficient functions for the quantities mentioned above Eq. (|L~l|), with the ex- 
ception of (x) in SIA, have been obtained in Refs. [2-7] (the latter two new articles deal with 
the only theoretically relevant Higgs-exchange structure function F§ in the heavy-top limit). The 
NNLO spacelike (S) splitting functions in Eq. (Of) are fully known from Refs. [8], while for the 
timelike (T) case only the diagonal quantities Pqq]g g (x) have been derived so far [9]. At N 3 LO only 
the coefficient functions for the structure functions fi.2,3,0 are available at this point [6, 10]. 

2. In ( 1 — x) contributions to the splitting functions and coefficient functions 



From order a 2 the quark coefficient functions in Eq. (1.2) and quark-quark splitting functions in 



Eq. (1.3) need to be decomposed into (large-* dominant) non-singlet and (suppressed) pure-singlet 



contributions. The non-singlet splitting functions receive an only single-logarithmic (SL) higher- 
order enhancement, and that only in terms relatively suppressed by (1— x) k ~ 2 [8,9, 11, 12] 

P$(x) = Az+^l-x); 1 +fl /+1 5(l-jc) +C/+i Iii(I-jc) + ^((1-x)*-V(1-jc)) . (2.1) 

Also the Cp = part of the gluon-gluon splitting functions is of this form [6]. The corresponding 
(pure-) singlet splitting functions include double-logarithmic (DL) contributions 
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p is'Ll l p m- terms up to (1-x) hr^l-x) / ln 2/ (l-x) . (2.2) 



' Ps,gglc F 



The non-singlet coefficient functions for the structure functions F\ ,2,3, the (transverse, angle- 
integrated and asymmetric) fragmentation functions Fjj,a and the quark-antiquark annihilation DY 
cross section F DY , on the other hand, show a DL enhancement already (but not only) at the (1 — x)^ 1 
plus-distribution level, i.e, 



c, 



(/) 



o.ns • 



terms up to (1-x) 1 \a 2l ~ 1 (l-x) . (2.3) 



The highest (1— x)~ l logarithms are resummed by the threshold exponentiation [13], with the ex- 
ponents now known to next-to-next-to-next-to-leading logarithmic (N 3 LL) accuracy (up to the nu- 
merically irrelevant four-loop cusp anomalous dimension^ in Eq. (p|)), see Refs. [14-16]. The 
leading contributions for the corresponding longitudinal DIS and SIA coefficient functions are 



down by a factor (1— x) and one power of In (1 — x) w.r.t. Eq. (2.3). Despite a recently renewed in- 
terest in such terms which behave as A -1 ln k N in Mellin space, see, e.g., Refs. [17], corresponding 
resummations have not been derived so far for these contributions to the coefficient functions. 

In the flavour-singlet sector we will confine ourselves to the DIS cases of F2 and F$ [6, 7] 
in the present contribution. The gluon coefficient function for F$ is of the form ( |2.3[ ), while the 
'off-diagonal' (see section 5) quantities are also double-logarithmic but suppressed by (1— x), 

c i\/^q : termsu P to ln 2 ' _1 (l-x). (2.4) 

3. Non-singlet physical kernels and coefficient-function predictions 

We now switch to moment space (and often suppress the Mellin variable A), which considerably 
simplifies the following calculations by turning the convolutions in Eqs. ( |ET| ) and ( 13) into simple 



products. The resulting manipulations of harmonic sums [18] and harmonic polylogarithms [19] 
have been mostly carried out in FORM 3 and TFORM [20]. 

The non-singlet physical evolution kernels K a for the DIS and SIA cases are constructed by 
^ d T2 ^ = -Jrkv + C ° P « = T 1 + C ° P ) C a lp a 



d\nQ 2 ding 2 v d\nQ 2 ' ' ' </« 



Pa + Pith) < ^ L ) p a = K a F a = £ /=0 a l + l K a jF a (3.1) 

for /I 2 = Q 2 (the additional terms for \i 2 ^ Q 2 can be readily reconstructed), where j3(a s ) is the 
usual beta function of QCD, /3 (a s ) = — a 2 j8o — a s 3 j8i — . . . with j3o = 11/3 Ca — 1/3 n* etc, and n* 
is the number of effectively massless flavours. For a^L Eq. ([Oh leads to the expansion 



K a = a s P afi + £ l=1 a' 5 +l (P a j - £^ p k c a ^ k ) (3.2) 

with 

Ca,\ = Q,l > C<3,3 = 3c a 3 — 3c a .2C ai i +C t 3 j 

c a .2 = 2c a;2 -c 2 i , CaA = 4c aA - 4c a j c a> \ - 2c a2 +4c 0)2 c al - c al , ... . (3.3) 

to 



The structure for the DY case is the same except for P a n = P a ns — > 2P a n in Eq. (3.2). 
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The threshold resummation of these coefficient functions [13], again for a ^ L, is given by 

C a (N,a 5 ) = g (a s ) cxp{L gl (a 5 L)+g 2 (a s L) + ...} + ff(l/N) (3.4) 



with L = InN. Due to the logarithmic derivative in the second line of Eq. (3. 1 ), the exponentiation 



(3.4) guarantees a single-logarithmic large-TV/Targe-x enhancement of the physical kernels [21], 



£ fl (Af,a s ) = -£ /=1 Aia l s L + P(a s ) J- {Lg l (a s L)+g 2 (a s L) + ...} + ... . (3.5) 

We are now ready to present the first crucial observation: all considered non-singlet kernels 
K a (including a = L in DIS and SIA) are single-log enhanced to all orders in N~ l or (1— x) [12,22]. 
Switching back to x-space, the universal a ^ L leading-logarithmic terms in DIS (upper sign) and 
SIA (lower sign) to N 3 LO read, with p qq (x) = (1 — x)~ ! — 1 — x, 

K a,o( x ) = 2C>p qq (x) + 3Cy<5(l-x) 

K al (x) = ln(l-x)p qq (x)[-2C Fj 8 T8C|lnx] 

K a2 (x) = ln 2 (l-x)/? qq (x) [2C Fj 8 2 ± 12C|j8 lnx+16C|ln 2 x] (3.6) 

K a3 (x) = ln 3 (l-x)/? qq (x) [-2C F /3 3 T44/3C| j 8 2 lnx-32C| j 8 ln 2 x + ^ 3 C^ln 3 xj , 



where % P is the unknown four-loop SL coefficient in Eq. (2.1). For DIS the N 3 LO relation is based 



on Refs. [10], while for SIA we have used incomplete but sufficient analytic-continuation results 



presented in Ref. [12] where also the DY relations analogous to Eqs. (3.6), known only to NNLO, 
can be found. The first terms on the right-hand-sides include the leading large-n^ terms for which 



Eqs. (3.6) can be generalized to all orders in DIS, using the C2, n s results of Ref. [23]. 

It is now rather obvious to conjecture that the physical evolution kernels receive only SL 
contributions to all orders in a s at all powers of n*. This implies an exponentiation (see section 4) 



of the coefficient functions beyond the (1 — x) + l terms. The emergence of the resulting fourth-order 



predictions can be illustrated by recalling the last relation written out in Eq. ( |3.3[ ), 



CaA = 4 ^,4 -4c a;3 c a ,i - 2c 2 2 + 4c fl;2 c 2 i - c*, . (3.7) 



SL DL.new DL, known for DIS/SIA 

,7,6,5/ 



I.e., the In ' ' (1— x) DL fourth-order contributions for F\ 23 and Fjja in Eq. (L2) need to cancel 
the corresponding terms from the known lower-order coefficient functions at all orders in (1— x), 
and consequently can be predicted from those results. We do not have the space here to give an 
explicit example of these predictions and their numerical size, but refer the reader to Ref. [12]. 
Note, however, that our results explain an old observation [25] for the highest 1 /N logarithms. 

Due to the universality of the leading terms in Eqs. ( |3.6| ), also for c^j^ in DIS and SIA the 
coefficients of the three highest logarithms are predicted, by the respective differences K 2 — K\ and 
Kj — Kj. The agreement of these predictions with those obtained from the quite different kernels 
Kl [22] - the above differences are of the order (1 — x)°, while the leading large-x terms of Ki are of 
the form (1 — xK 1 - provides a quite non-trivial check of the above conjecture. On the other hand, 
only two logarithms can be predicted completely at this point at the third and all higher orders for 
the DY case [12], as the corresponding coefficient function is only known to order a 2 [5]. 



4 



Higher-order predictions from physical evolution kernels 



A. Vogt 



4. All-order exponentiation of the 1 /N non-singlet coefficient functions 

The subleading 1 /N contributions to the non-singlet coefficient functions for Fi.2,3, Ft ia and F^y 
can be cast in an all-order form analogous to (if unavoidably less compact than) Eq. (|3~4|), 

C a ~ C a \ N Q Lk = ^ ( [dj$L + ] As + K ( j L + d$] al + . . . ) exp {Lh, (a s L) + h 2 {a s L) + ...}. 

(4.1) 

The exponentiation functions are defined by the series hk(a s L) = Y,k=\ hu n {a s L) n with L = \nN. 
Their coefficients for DIS/SIA (given by the upper/lower sign in Eqs. ( |12| ) and (p~3|)) relative to 
the corresponding coefficients for the N°L k soft-gluon exponentiation are given by 

hk = Sik , h 22 = g 22 + Iflo 2 ± jPoC P - 18C| (4.2) 

hu = g 2l + ^ Q ±6C F , h 23 = £ 23 + ^A) 3 ±(^-§)i8 2 C f -y/3 C F 2 ±72C| . 



Note that only the Cfj3 ' and Cp j8 ' 1 terms of K u j in Eqs. (^6) are relevant at this order in 1 /N. 



^ is the corresponding subleading large-«y coefficient at the fourth order, the calculation of which 



should become possible in the not too distant future. Also the first term of h 3 in Eq. ( pklD is known 



for DIS and SIA but non-universal, as the effects of Fl set in at this point. See again Ref. [12] for 



these results as well as the prefactor coefficients in Eq. ( |4.1| ) and all corresponding DY results. 

The corresponding exponentiation for the longitudinal structure function and fragmentation 
function is given by [22] (see also Ref. [24]) 

(AO = N-\d ( L fa s + d { L fa 2 s + ...) exp {L/i La (a s L) +h L ^ 2 {a s L) + ...} + ^(A^ 2 ) , (4.3) 

where the following coefficients can be determined from the third-order result of Refs. [3, 12]: 

2 1 

^z.,11 = 2c f , h L l2 = -poC F , h L13 = -PqC f 

h L .2i = Po + 4 J e C F - C F + (4 - 4 (Q - 2C F ) 

h,22 = ^( j 8 / J21 +A 2 )-8(C A -2C F ) 2 (l-3C 2 + C 3 + C2 2 )- (4-4) 

Both these exponentiations have far less predictive power than their N°L k counterparts [13-16] 
where, e.g., unlike in Eqs. (4.2) and ( |4.4[ ), no other new coefficient enters g 22 besides the two-loop 



cusp anomalous dimension A 2 . A full NLL accuracy, i.e., a complete determination of the function 



h 2 {a s L) may be feasible for Eq. (11). On the other hand, the corresponding leading coefficients 
for h 3 [12] and the results for h L2 in Eq. ( |14| ) indicate a major, possibly insurmountable obstacle 
on the way to full NNLL and NLL accuracy for the quantities F a (a ^ L) and Fl, respectively. 

5. The singlet evolution of F2 and and splitting-function predictions 

DIS via the exchange of a scalar (p directly coupling only to gluons (like the Higgs boson in the 
heavy-top limit [26]), is an ideal complement to the standard structure function F 2 . The evolution 



kernels for the resulting system of observables are as in the first line of Eq. (3.1), but with 



ft)' c = ( C r' q r' 8 )> k =( K k 22 k 2 ^ i5A} 
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and the splitting-functions matrix Py. This system has first been discussed at NLO in Ref. [1] (it 
may also be interesting to study other systems such as (Fz, F£) [27] and corresponding SIA cases). 
Instead of the second line of Eq. (fTl]), we now have (with [C,P] denoting the matrix commutator) 



dF 
d\aQ 2 



d\nC 




+ [C,P]C- l +P)F = KF . 



DL (singlet only) 



(5.2) 



As far as they are completely known now, i.e., at NLO and NNLO, also the matrix entries of 

(n) 

K show an only single-logarithmic enhancement at all powers of (1— x), K J ~ ln"(l— x) + 

(n) 

Moreover, the leading-log contributions to ^22/00 316 tne same as m tne non-singlet quark-case 
and the very closely related Cp = gluon case [6]. Conjecturing that this behaviour holds also at 
N 3 LO, the highest three logarithms of the unknown four-loop splitting functions, 



In 



6,5,4, 



1- 



of P 



(3) 



qg,gq 



and ln 5 ' 4 ' 3 (l-x) of P {3) . 

v ' Ps.gglcp 



(5.3) 



can be predicted from the known [6, 10] three-loop coefficient functions for F2 and F§ at all orders 
in (1 — jc). For example, the leading (1— x) part of the N 3 LO gluon-quark splitting function reads 



o (3) (x) = ln 6 (l-x) -0 + ln 5 (l-x) 
+ ln 4 (l-x 



22 ^ 14 ^ ^"2 2 

— Op/lf - — Cj^CFtlf - — CfipTlf 



/ 293 
LV127" 



80.^3 ^/4477 „ r \ r2r 



81 



9 

116 2 2 17 2 4 3 

"777" Qf n / + 777 Qf^"/ — 777 Qf w , 
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51 



162 

X 

8 J c w 



■ CipCpUt 



+ ^(ln 3 (l-x)) (5.4) 



with = Ca — Cf. The vanishing of the leading In 6 (1— jc) term is due to an accidental cancella- 
tion of positive and negative contributions to its coefficient. On the other hand, the colour factors of 



the DL terms in Eq. (5.4) follow the same pattern as the corresponding lower-order contributions: 
all DL terms vanish for Ca = Cp (part of the supersymmetric limit), with the leading terms of Pqg 
being of the form n^C^p, the next-to-leading logarithms n^{Cp ,nAC^ etc. Rather non-trivially, 
this pattern is predicted to hold for all four-loop singlet splitting functions at all orders in (1— x) [6]. 

Unlike in the non-singlet case, there is no direct all-order generalization here, as the cancel- 



lation of the DL contributions in Eq. ( p.2| ) involves the corresponding terms ( |2.2[ ) and ( |2.4[ ) of the 
N"LO splitting functions and coefficient functions which are both unknown at n > 4. One may 
try a simultaneous extraction of at least the leading logarithms of both quantities, but it turns out 
that the only single-logarithmic enhancement of the physical kernel does not quite provide enough 
constraints, even if the colour structure of the previous paragraph is assumed in addition. 
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